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ABSTRACT 

We propose that particles during the second phase of solar 
flares are accelerated by stochastic resonant scattering 
off hydromagnet ic waves and first-order Fermi acceleration 
in shock waves generated in the impulsive phase of the 
flare. Our solutions allow arbitrary power law momentum 
dependences of the momentum diffusion coefficient as well 
as the momentum loss time. The acceleration time scale to 
a characteristic energy 'dOOkeV for protons can be as short 
as 5s. The resulting electron spectra show a characteristic 
double power law with a transition around 200 keV and are 
correlated to the proton spectra evaluated under equal 
boundary conditions, indicating that electrons and protons 
are accelerated by the same mechanism. The correlation 
between the different spectral indices in the electron 
double power law and between electron and proton spectra 
are governed by the ratio of first-to second-order accel- 
eration and therefore allow a determination of the Alfven 
Mach number of the shock wave. 


1 . Introduction. We propose that the combination of first- and second- 
order Fermi acceleration is the mechanism responsible for second-phase 
acceleration in solar flares. This model has a number of distinct advan- 
tages over previous models, in that it gives a natural explanation for 
the various observed time delays between first and second phase as well 
as for the dependence of the particle spectra to the strength of the flare 
and the correlation between electron and proton spectral indices. 


2. Theory. Parker and Tidman (1958) have pointed out that Fermi accel- 
eration is an intrinsic property of any sufficiently agitated plasma of 
energetic particles. The effects of second-order Fermi acceleration due 
to irregularly moving magnetized fluid elements, first-order Fermi ac- 
celeration off strong shocks as well as loss and escape processes can be 
incorporated into a transport equation in phase space (Ramaty 1979) 
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where p is the particle momentum, N(p) = 4up 2 f(p,t) the number of par- 
tides per unit momentum and unit volume, and Q(p,t) represents sources 
and sinks of particles. Using the concept of the age distribution 
(Schlickeiser and Lerche 1985) the effects of spatial diffusion, convec- 
tion and catastrophic losses have been combined in a momentum-dependent 
escape time T(p) = T Q p” b . 


In a plasma with a strong MHD turbulence both the momentum diffusion 
coefficient D(p) and the spatial diffusion coefficient along the mean 
magnetic field K„(p) are governed by the turbulence simultaneously. By 
using quasi-linear theory K„ (p) can be related rigorously to the spectrum 
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of the magnetic field fluctuations (Jokipii 1977). Assuming a power law 
spectrum for the magnetic irregularities W(k) = W Q k - 9 we may write 

Kn(p) = <5 v(p) p2 - q = Kpl (2) 

(6,K = const.) and 

°( p ) = a 2 k7 P 2 

(V^: Alfven speed, a 2 : constant). As gain process, we consider quasi- 
continuous momentum gain by acceleration at a shock wave moving through 
the plasma with speed Vg. The momentum gain by first-order Fermi accel- 
eration at a single shock wave has been determined by Drury (1983) 
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(a x : constant). Mobius et al. (1982) pointed out that for conditions in 
the flare site the rate of systematic acceleration exceeds the rate of 
momentum loss. Thus we assume that pg >> p^ at least at all momenta of 
interest. Using (2), (3) and (4) in equation (1) and considering the 
steady-state case (3f/8t = 0) yields 
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(a 1 =a 1 Vg/K, a 2 =a 2 V^/K, a = a x /a 2 , A = 1/a 2 T ). We assume that the in- 
jection takes place at some momentum Q(p) = q Q o(p-p Q ) and we find for the 
steady-state particle number density N(p) = 4 it p 2 f(p) (see Droge and 
Schlickeiser (1985) for details): 
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V = I (3+a-n)/(u+b) 
law distribution 


In the limit n+b •> 0 solution (6) approaches a power 
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Equations (6) and (7) are generalizations of the solutions of Ramaty 
(1979; b = 0, a = 0, r| — 0, 1 ) , Barbosa (1979; a = 0) , Pikel'ner and Tsytovich 
(1976; a = 0, b = 0) . The solutions (6) and (7) allow arbitrary momentum 
power law dependences of the spatial diffusion coefficient (n) and the 
escape time (b). The parameter a = (a 1 /a 2 ) (Vg/V^) = (a x /a^) is of the 
order of the square of the shock's Mach number and describes the ratio of 
first- to second-order Fermi acceleration. In the case of no Fermi ac- 
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Fig. 1: Normalized steady-state 
electron number density. Parti- 
cles are injected at E o =50keV. 
At E t = 200 keV a transition 
from solution (6) to (7) occurs. 
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Fig. 2: Normalized steady-state 
proton number density. Particles 
are injected at E Q = 10 keV. 


celeration at shock waves (a * 0) (6) and (7) reduce to the solutions of 
previous models as quoted above. 


3. Discussion. The energy spectra of electrons accelerated in large 
flares exhibit a characteristic double power law with a break around E T s 
200 keV (Lin et al. 1982). From this we conclude that for nonrelativistic 
kinetic energies (E < Ex) the Bessel function solution (6) holds, which 
for small arguments is approximately a power law, whereas for E > Ej the 
particles are relativistic (p=mc=E = 200 keV for electrons) and the 
distribution function approaches the power law solution (7) (cf. Fig. 1). 
We take spatial diffusion along the mean magnetic field as the relevant 
escape process. The escape time then is T(p) = L 2 /K„ = (L 2 /k) p - W, where L 
is the length scale of the system and n = b = 0 in the relativistic and 
n =b = 1 in the nonrelativistic case. Thus we may consider the evolution 
of the particle spectra in the hard-sphere approximation (cf. Mobius et 
al. 1982), which gives us 

K„(p) = y v(p) it (8) 


v(p) is the velocity and it the momentum independent mean free path of the 
particles. For comparison with data we transform (6) and (7) into energy 
space where N(E) =N(p) (1/v) and E = (p 2 c 2 +m 2 c‘*) 1/2 -me 2 is the particle 
kinetic energy. An "effective power law" exponent can be calculated (E >E 0 ) 
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Fig. 3: Correlation between elec- 
tron spectral exponent below 200 
keV Yt and exponent above 200 keV 
Y_j as a function of a 
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Fig. 4: Correlation between Yj 
and proton spectral exponent 
above 10 MeV as a function of a 


x = £ c /LVa is a free parameter which may have different values for elec- 
trons and protons and for different flares. Equations (9) and (10) allow 
us to determine the correlation between the mean spectral index below 200 
keV Yl and above 200 keV Yu for electrons (Fig. 3) as well as the cor- 
relation between Yi and the proton spectral exponent above 10 MeV Y p 
(Fig. 4) as a function of a. The curves obtained from (9) and (10) are 
in excellent agreement with the measurements of Lin et al. (1982) , indi 
eating that the strongest shock waves have an Alfven Mach number M^ 'v 4, 
corresponding to a n* 30 (see Droge and Schlickeiser (1985) for details). 

4. Conclusions. Combining first- and second-order Fermi mechanism in 
solar flare second— phase acceleration successfully explains the observed 
ion and electron energy spectra. The model correctly accounts for the 
sometimes very short delay times, and reproduces quantitatively the cor- 
relations of nonrelativistic with relativistic electron spectral indices, 
and the correlation of nonrelativistic electron with proton indices. 
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